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ABSTRACT 

Modified Newtonian dynamics (MOND) proposed by Milgrom provides a paradigm 
alternative to dark matter (DM) that has been successful in fitting and predicting the 
rich phenomenology of rotating disc galaxies. There have also been attempts to test 
MOND in dispersion-supported spheroidal early-type galaxies, but it remains unclear 
whether MOND can fit the various empirical properties of early-type galaxies for the 
whole ranges of mass and radius. As a way of rigorously testing MOND in elliptical 
galaxies we calculate the MOND-predicted velocity dispersion profiles (VDPs) in the 
inner regions of ^ 2000 nearly round SDSS elliptical galaxies under a variety of as¬ 
sumptions on velocity dispersion (VD) anisotropy, and then compare the predicted 
distribution of VDP slopes with the observed distribution in 11 ATLAS^^ galaxies 
selected with essentially the same criteria. We find that the MOND model parame- 
terised with an interpolating function that works well for rotating galaxies can also 
reproduce the observed distribution of VDP slopes based only on the observed stellar 
mass distribution without DM or any other galaxy-to-galaxy varying factor. This is 
remarkable in view that Newtonian dynamics with DM requires a specific amount 
and/or profile of DM for each galaxy in order to reproduce the observed distribution 
of VDP slopes. When we analyse non-round galaxy samples using the MOND-based 
spherical Jeans equation, we do not find any systematic difference in the mean prop¬ 
erty of the VDP slope distribution compared with the nearly round sample. However, 
in line with previous studies of MOND through individual analyses of elliptical galax¬ 
ies, varying MOND interpolating function or VD anisotropy can lead to systematic 
change in the VDP slope distribution, indicating that a statistical analysis of VDPs 
can be used to constrain specific MOND models with an accurate measurement of 
VDP slopes or a prior constraint on VD anisotropy. 

Key words: dark matter - galaxies: elliptical and lenticular, cD - galaxies: kinematics 
and dynamics - galaxies: structure - gravitation 


1 INTRODUCTION 


Mass discrepancy (or missing mass) problems in galax¬ 
ies, galaxy clusters and the Universe have usually been 
attributed to dark matter (DM) assuming that standard 
Newton-Einstein gra vity is valid in those dynamical systems 
(for r e views, see, e.g., TrimbldllQSTI : iBertone. Hooper fc SilkI 
I2OO5I : ISandersI 1201 4al l. The currently popular Lambda 
cold dark matter (ACDM) cosmological paradigm of the 
Universe provides a successful phenomenological model 
of the large-scale structure, the cosmological microwave 
background radiation (CMBR) anisotropy power spec¬ 
trum, and the expansion history of the Universe (e.g. 


iHinshaw et al.l l2013l : lAde et ^ l2014l l. Theories of galaxy 
formation and evolution based on the ACDM paradigm 
have also been making progress with the goal of explain¬ 
ing the rich phenomenology of galaxie s (for a review see, 
e.g., IMo. van den Bosch fc Whitg|2oi^ . 

However, in the realm (where dynamics is non-linear) of 
galaxies the ACDM paradigm appears to face serious chal¬ 
lenges. Rotation curves of disc galaxies start to deviate from 
Newtonian expectation all at the same critical acceleration 
scale ao ~ 1.2 x 10“^°m s“^ regardless of galaxy size, lumi¬ 
nosity, surface brightness or any other property (iMcGauehl 
I2OO4I : iFamaev fc McGauehl [20121 . The velocity at the flat 


© 2015 RAS 

























2 Chae & Gong 


part of the rotation curve is well correlated with the total 
baryonic mass of the galaxy (consistent with zero intrinsic 
scatter) spanning five orders-of-magnitude in baryonic mass 
from giant spiral galaxies to low surface brightness galaxies; 
this corr elation is now known as the baryonic Tully-Fishe r 
relation llMcGaughll2005l . I 2 OIII : l^maev fc McGaugljl2012h . 
Furthermore, to every feature in an individual rotation curve 
there corresponds a feature in baryonic mass density pro¬ 
file even in galaxies whose dynamics should be dominated 
by DM if the ACDM par adigm is assumed (ISancisil l2004l : 
iFamaev fc McGaughll201^ 1. 

The apparent existence of the critical acceleration no 
and the intimate connection between the galaxy rotation 
curve and the baryonic mass distribution seems unnatural 
in the ACDM paradigm. Other outstanding challenges of 
the ACDM paradigm in galaxies such as the number and 
distribution of DM subhaloes and the DM central cusp/core 
density profile are also we l l documented in the literature 
jFamaey^jgMdhGnugd iKroupa, Pawlowski fc MilgromI 


l2012l:lMcGaughll2014li. 


Milgrom’s modified Newtonian dynamics (MOND; 

I MilgromI 1 1983h offers an alternative paradigm that can by¬ 
pass the above challenges of rotating galaxies faced by the 
ACDM paradigm. The MOND paradigm posits that New¬ 
ton’s law of gravity breaks down at extremely low accel¬ 
erations below the critical acceleration oq. Then, the phe¬ 
nomenological existence of oq, the baryonic Tully-Fisher re- 
lation, and so on can be u nified in a single MOND law (see 
iFamaev fc McGaughl20i^ for an extensive review). Whether 
the empirical MOND law can be ultimately explained by a 
fundamental law of modified gravity or somehow by revised 
theories of galaxy formation under the ACDM paradigm, it 
is clearly worthwhile as a Kepler-like law of galactic kine¬ 
matics. If the remarkable success of the MOND paradigm in 
rotating galaxies implies an underlying law, it must also ex¬ 
plain the empirical properties of spheroidal early-type (ellip¬ 
tical and lenticular) galaxies (ETGs) that are different from 
disc galaxies in dynamics and structure. In this respect it is 
quite interesting to test MOND in ETGs. 

However, testing MOND (or DM) in ETGs is far more 
challenging because they lack neutral hydrogens (the cru¬ 
cial dynamical tracer for disc galaxies) and orbits of stars 
(and other tracers) are complex (and/or unknown). Never¬ 
theless, significant efforts have been made to test MOND 
(and DM) in ETGs using dynamical tracers such as stars, 
planetary nebulae (PNe), globular clusters (GGs), hot X- 
ray gases, satellite galaxies, and gravitational leasing. De¬ 
spite significant efforts, it is not yet clear whether MOND 
can provide a successful phenomenological model for the dy¬ 
namics of ETGs. Sometimes, analyses of the same galaxies 
often reached different conclusions. 

Based on the kinematics of PNe (out to ~ 5Re 
whe re Re is the effectiv e radi us) of three elliptical galax¬ 
ies, iMilgrom fc Sand^ (l2003l ~l argued that MOND could 
explain the observed kinematics well an d no DM would be 
needed. However, for the same galaxies IPekel et al.l (I 2 OO 5 I I 
argued that DM would be consistent with the data once 
radially varying velocity dispersion (VD) anisotropies were 
allowed. Independ ent analyse s of P Ne kinematics of two 
other ellipticals bv lTiret et al.l (120071) and lSamurovid (I 2 OIOI I 
showed that MOND could explain the kinematics without 
DM. Based on the velocity dispersion profiles (VDPs) of 


stacked SDSS satellite galaxi es around two n a rrow lumi¬ 
nosity range s of re d galaxies iKlvpin fc Pradal (120091) and 
lAngus et al.l (I2OO8I ) reached opposite conclusions respec¬ 
tively against and for MOND because of their different as- 
sumption s on the profile s of density and VD anisotropy of 
satellites. iMilgronil (l2012f) tested MOND using the mass pro¬ 
files of two X-ray bright elliptical galaxies deduced assum¬ 
ing hydrostatic equilibrium and found that MOND success¬ 
fully reproduced the profiles over a wide acceleration range 
from > lOtto (near the galactic centre) down to ~ O.lao 
(over 100-200 kpc). Analyses of (the masses of) strong lens 
ETGs also pr oduced confusing res ults, some arguing against 
MO ND (e.g. Fen:erasefL_aI, 2012ll wh ile others for MOND 
(e.g. IChiu et al.l 20 111 : Sandersl 2014lJ) . 

Based on t he kinematics of GGs (out to several Re) 
several studies (lSamurqyic_&^_Cirk^^ 200^ RIcMer et al.l 


I2OO8I : ISchuberth et al.ll2012l : ISamurovid 12012 . 120l4) tested 
MOND in a han dful of elliptical galaxies. In particular, 
ISamurovij (|2014|) carried out a systematic analysis of ten 
elliptical galaxies covering a broad range of mass. These 
studies of GG kinematics employed spherical Jeans analy¬ 
ses assuming isotropy or rel atively simple ani sotropies. Ac¬ 
cording to these studies, as ISamuro\^ (I20l4) emphasized, 
MOND alone without an additional DM component could 
not provide a successful fit for several massive slow rotators 
in their outer regions (beyond ~ 3Re), although MOND was 
generally successful for less massive f ast rotators (where th e 
slow/fast dichotomy is in the sense bv lEmsellem et ^l2007l ) . 
iRichtler et al.l (l2008l) also advocated the need of DM under 
MOND in the central massive galaxy of the Fornax cluster. 
It is unclear whether these results imply the breakdown of 
MOND for a class of elliptical galaxies. The main limits of 
these studies are small sample sizes and uncertainties in VD 
anisotropy. 

A rigorous test of MOND in ETGs would require (1) a 
statistically representative sample of ETGs and (2) an em¬ 
pirical probe of mass profile over a wide acceleration range. 
Here, we propose a test of MOND in elliptical galaxies sat¬ 
isfying the first requirement to a large extent and the sec¬ 
ond requirement only partially. We calculate MOND-based 
VDPs of ~ 2000 nearly round galaxies from Chae, Bernardi 
& Kravtsov (2014) for 0.1J7e < r < Re using the observed 
stellar mass distributions under a variety of possibilities of 
VD anisotropy. We calculate the slopes of the VDPs for 
0.1i?e < r < Re and compare the distribution of the slopes 
with the measured distributio n in a similarly selected sample 
of 11 nearly round ATLAS®® (iGappellari et al!ll2011^ galax¬ 
ies. Hence our analysis is statistical in nature. We describe 
the data and the method of analysis in section 2 and the 
results in section 3. We discuss implications of our results 
for MOND and give our conclusions in section 4. 


2 DATA AND METHOD OF ANALYSIS 

We consider only elliptical galaxies that are nearly round so 
that they can be analysed based on the spherical symmetry 
assumption with minimal error. A statistically representa¬ 
tive sample of nearly round galaxies can b e drawn from the 
immense Sloan Digital Sky Survey (SDSS: IVork et ^l2000f) 
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data baseQ We use a sample of SDSS g alaxies defin e d and 
analyse d under the ACDM p aradigm bv IChae et alJ ll2014l i 
(see also lBernardi et '^l201(j ) that contains pho tometric pa¬ 
ramet ers re cently invest i gated and measured bv iMeert et alJ 
ll2013li and iMeert et alJ (I2015l i based on data release (DR) 
7 data0 The sample contains 2054 elliptical galaxies with a 
mean redshift of (z) ~ 0.12 spanning two orders of magni¬ 
tude in stellar mass that are nearly round (projected minor- 
to-major axis ratio b/a > 0.85) and do not possess a mea¬ 
surable disc. 

Because each galaxy has no measurable disc by selec¬ 
tion, its light distribution on the sky can be well-described 
by a Sersic (1968) profile 


where p{r) is the volume luminosity density at radius r, (Jr{r) 
is the radial stellar VD, and /3(r) is the VD anisotropy given 
by /3(r) = 1 — ut{r)/afir) where o't(r) is the tangential VD 
in spherical coordinates. In equation ([2]) g{r) is the MON- 
Dian (‘effective’ or ‘real’) gravitational acceleration which 
differs from the Newtonian acceleration gN{r) = GMh{r)/r^ 
for the baryonic (stellar here) mass Mh{r) within r. The ac¬ 
celeration g approaches gN for g/ao 1 but tends to y/gWoo 
as g/ao 0 (hence g{r) due to a point mass makes a tran¬ 
sition from the Newtonian inverse square law g{r) oc r~^ 
for g/ao 1:^ 1 to g{r) oc r~^ for g/ao <C 1). The transi¬ 
tion between the Newtonian and the MONDian regimes is 
parameterised by 


E(R) oc exp 


t>nser 



( 1 ) 


where nser is referre d to as Sersic index and = 2nser — 
1/3-1-0.009876/nser llPrugniel &: Simi^ll997^ . Deprojecting 
equation © we obtain a volume luminosity density profile 
p{r) and then a stellar mass density profile through a stel¬ 
lar initial mass function (IMF). (For our selected galaxies 
gas is ignorable particularly because they do not possess 
a measurable disc.) Each galaxy has also a measured VD 
Uap within the SDSS aperture of Rap = 1.5 arcsec radius 
which corresponds to a physical radius between O.lRe and 
Re with a mean of (Rap/Re) ~ 0.5. Note that Uap denotes a 
luminosity-weighted average of line-of-sight velocity disper¬ 
sions (LOSVDs) of all stars within the aperture. 

Under the standard Newtonian dynamics the observed 
VD CTap may be reproduced by adjusting the stellar mass- 
to-light (M*/L) ratio or equivalently the stellar IMF for a 
flexible assumption of DM mass profile (even including the 
case of no DM) given the empirical intrinsic galaxy-to-galaxy 
scatter of IMF. This means that Uap (an average quantity 
within a region) itself may not be a sensitive probe of grav¬ 
ity or mass distribution. However, the VDP is a powerful 
probe. The observed VDP of an individual galaxy over a ra¬ 
dial range can be used to infer the tota l mass profile for tha t 
range under Newtonian dynamics ('e.g. lThomas et al]l2007li . 
For a large sample of galaxies that do not have i ndivi dually 
measured VDPs one can, as done in IChae et al.l ll2014l ~). pre¬ 
dict VDPs using galaxy models satisfying the available ob¬ 
servational constraints, and then compare the statistical dis¬ 
tribution of VDPs with an available empirical distribution. 
Such a comparison of VDP slope distributions clearly shows 
that galaxies must be embedded in DM halos to match the 
empirical distribution if Newtonian d ynamics is assume d to 
be valid at all acceleration scale fsee lChae et al.ll2014l and 
below). 

Here we carry out a similar statistical analysis of VDP 
slopes under MOND. Our analysis is b ased on the spheri¬ 
cal Jeans equation giv en in MOND (e.g. lAngus et al.ll200a : 
iKlvoin fc Pradall2009l i by 


d-lp{r)c^fir)] 

dr 


+ 2^[p{r)afir)] 
r 


-p{r)g{r), 


( 2 ) 


^ http:// www.sdss.org/ 

^ The iMeert et alJ <l2015l^ photometric measure¬ 
ments for full DR7 data can be downloaded at 
http: //shalaowai.physics.upenn.edu/'^ameert/fit_catalog/download/ 


gN = p{g/ao)g (3) 

where interpolating function p,{x) satisfies p,{x) —^ 1 for x 
1 and p.{x) —>■ X for a; ^ 1. The transition can also be written 
as the inverted relation 


g = v{g^/ao)gn, 


(4) 


where iy{y) satisfies u{y) 1 for y ^ 1 and I'iy) y~^ ^ 
for j/ <C 1. 

We consider a class of interpolating functions given by 


Unix) 


X 

(1 -|-a:")V" 


(5) 


with _ the _ correspo nding inverted function 

llMilgrom fc Sand^l2008l i 


vniy) = 


1 -f (1 -f 

2 



( 6 ) 


where the case n = 2 is traditionally known as the ‘stan- 
dar d’ function while t he ‘sim ple’ case n = 1 introduced 
by iFamaev fc Binn^ ll2005h has turned out performing 
well i n various recent studies (e.g. ISa^ers fc Noordermeeii 
I 2 OO 7 I : lAngus et al.l l2008l : iMilgromI l2012l i. We also consider 
the interpolating function implied by Bek enstein’s relativis¬ 
tic theory of modified gravity given by llBekensteirJ l2004l : 
IZhao k, Famaevll2006li 


, . VI + 4a: — 1 
peek (a;) = , —- 

VI + 4a: -I-1 

with the corresponding inverted function 

i'Bek(t/) = 1-l-y”^^^. 


(7) 

( 8 ) 


For MONDian dynamical analyses of a galaxy any ex¬ 
ternal field in which the galaxy is embedded should, in prin¬ 
ciple, be t aken into account. A MONDian external field ef¬ 
fect (EFE: lMilgr ornl ll983l : iFamaey. Bruneton fc Zha3l2007l : 


iRichtler et al. 2011 1 may particularly matter for those in 
the central regions of galaxy clusters and those having close 
neighbors. Most galaxies in our SDSS sample are not in clus¬ 
ter centres and have light distributions “uncontaminated” 
by neighbors. Also, an evaluation of EFE s for elliptical 
galaxi es using the Virgo and Coma clusters bv iRichtler et al] 
ll201lll indicates that the expected effects are small. In this 
work we do not consider including EFEs for our selected 
galaxies. 

For a given baryonic mass profile Mb(r) with the corre¬ 
sponding luminosity density p(r) the solution of equation ([2ll 
for the radi al stellar VD a^jr ) can be given following ap¬ 
pendix B of IChae et al.l (l2012l i as 
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(Jr (r) = G 


f 


fGA^\ 

u}{r) p{r) \ao ) ’ 


(9) 


where w(r) = exp [J’^(2/?(r')/r')dr'] for an anisotropy 
and vipy) is an inverted interpolating function. The LOSVD 
of st ars at projected radius R on the sky is then given by 
(e.g. iBinnev fc Mamonlll982h 


O'los(-R) = 


E(i?) 


roo 

Jr^ 


p{r)a^{r) 


p2 


dr^ 


Vr2 - i?2 


( 10 ) 


where E(i?) is given by equation ([T]). The luminosity 
weighted LOSVD within projected radius R is given as 


a{R) = {mos)iR) = 


ff E(R')aio4R')R'dR' 

ffS(R')R'dR' 


( 11 ) 


and thus the aperture VD is = o'(R — i?a,p)- 

Three factors can be involved in the connection between 
the observationally derived luminosity density p(r') and the 
observed VD (Japi the stellar mass-to-light ratio Mi,/L, the 
VD anisotropy profile /3(r) and the (inverte d) interpolating 
function i/loxlr'l/anL iBernardi et al.l (l201Cll l provide stellar 
masses (or stellar mass- to-light ratios) of SDSS ETGs based 
on the IChabrierl (l2003l ~l IMF. Following the literature (e.g. 


IChae et al.l 2014 ) we define a stellar mass-to-light ratio or 
IMF mismatch parameter 


5m = logio(M*/M*^‘') = logio 


(M*/L) 


L(M*/L)ChJ 


( 12 ) 


between the unknown stellar mass M* (or Mi,/L) and 
the hducial st e llar in ass [or deriv ed by 

iBernardi et al.l ll2010ll based on the Ichabrierl ll2003ll IMF. 
All the uncertainties related to photometric measurement, 
mass-to-light ratio, and/or IMF are then absorbed into this 
single parameter Sm- Since our goal is to test MOND i.e. 
its interpolating functions, we need empirical inputs or as- 
sumpti ons for Sm and /3(r ) 0 Here we take an approach sim¬ 
ilar to IChae et al.l ll20l4l . The empirical input for Sm is 
provided b y the literature results o n VD-dependent stellar 
IMF (e.g. ^ppellari et al.ll2013bl: IConroy fc v an DokkumI 


l2012l : lTortora. Romanowskv fc Napolitan oli2013l: see below). 

For /3(r) we c onsider both const a nt /3 and do uble Osipkov- 
Merritt-type dChae et al.l l2012l : IChael I20l4) /IdOM(r) = 
/3i/(l-|-rf/r^)-|-/32/(l-l-ri/r^), where Pi+I32 = Pir ->• oo) = 
Poo- For constant P each value is drawn from a Gaussian 
distribution. For 5dOM(r) the radial mean for 0 ^ r (g: Re, 
referred to as /3elj and the value at infinity {Poo) are drawn 
from a Gaussian distribution and ri and r 2 are randomly as¬ 
signed satisfying 0 < ri/Re < r^jRo < 1. The adopted func¬ 
tion /IdOM(r) is intended to encompass (not exhaustively) 
possible behaviours of anisotropy. 

For input distributions of Sm and P{r) (constant or 
double Osipkov-Merritt-type) we search randomly for a 
pair reproducing Uap within a typical measurement error of 
0.04 dex for the observationally derived p{r) of each galaxy. 
For most galaxies a pair is easily found. For a few per cent of 


® If a VDP is measured for a radial range of a galaxy, then one 
can in principle simultaneously solve for Sm and P{r) for that 
range for the measured p(r) and a(R). 

The relatio n between (de, Poo ) and {Pi, P 2 ) can be found in 
appendix B of IChae et alj ll2014l i 


a{R)/ 



Figure 1. Profiles of luminosity-weighted LOSVD cr(R) (equa¬ 
tion 11) for 11 nearly round ATLAS^^ galaxies for which 
LOSVDs have been measured over 300 bins in the optical region 
(see the text). Red and blue profiles represent respectively slow 
and fast rotators. The abscissa is radius R on the sky normal¬ 
ized by the effective radius Re while the ordinate is normalized 
by (Je 2 = o-(R = Re/2). Small circles represent the mean of the 
displayed profiles. 


galaxies a successful pair could not be found after a signifi¬ 
cant number of trials. For each successful model we calculate 
a VDP o-{R) (equation [TT} as a function of radius R on the 
plane of the sky. We use a power-law approximation 


^(R) ^ { R y 

(Te2 \Re/2j 


(13) 


for O.lRe < R < Re, where (Je 2 = (j(R = Re/2) and the 
value of Tj for each model is determined through a least- 
square fit. The statistical distribution of rj for the galaxy 
models is then compared with an empirical distribution. 

Published values of rj for dozens of FTGs have a 
mean of {rj) = —0.06 ± 0.01 with a ga laxy-to-galaxy in- 
trinsic scatter of (j„ = 0.03 — 0.04 (e.g. iGhae et al.l 1201 j : 
ICappellari et ahl l2nn(t >. For a more reliable comparison 
with our n early round galaxies we select galaxies from the 
ATLAS®° (ICappellari et all I 2 OIII I public datcQ using the 
same criteria. Out of 260 ATLAS^^ FTGs 3 3 galaxies satisfy 
projec ted h/a > 0.85 in both photometric dCappellari et al.l 
l2013al l and kinematic (iFmsellem et al.ll2011^ distributions 
within both Re/2 and Re- Seventeen of them are elliptical 
and out of these fifteen have measured values of LOSVD over 
300 bins (pixels) on the plane of the optical region allowing 
reliable measurement of rj. Finally, excluding galaxies hav¬ 
ing the mean of bulge and to tal Sersic indices nser < 2.5 as 
for the SDSS ETG sample (IChae et al.1 12014^ we are left 
with 11 galaxies: they are NGC 3193, 3379, 4168, 4278, 
4283, 4374, 4458, 4552, 5173, 5638, and 5846. For each 
galaxy about 20 concentric rings are dehned within Re. 
LOSVDs of the bins in each r ing are averaged a nd weighted 
with the surface brightness llScott et al.ll20li^ of the ring 
and then these luminosity-weighted average LOSVDs give a 
(j(R) profile. The derived (r{R) profiles for the 11 galax- 


® http://www-astro.physics.ox.ac.uk/atlas3d/ 
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Figure 2. Dark matter fraction within the sphere of radius equal 
to Re (/dm) in SDSS and ATLAS®® nearly round galaxies as a 
function of stellar mass (M*), effective radius {Re) and Sersic 
index (nser)- Large dots represent 11 ATLAS®® galaxies. SDSS 
galaxies are represented by small dots with thick solid curves 
representing the median values, thick and thin dashed curves the 
68% and 95% limits respectively. See the text for details. 


ies are displayed in Fig. [T] For these galaxies the least- 
square fit power-law slopes for 0.1i?e < R < Re have a 
mean of (v) = —0.057 ± 0.011 with a standard deviation 
of 0.037 ± 0.007 (where the quoted errors have been es¬ 
timated from bootstrap resampling) con sistent with pub¬ 
lished results for ETG s of any ellipticity dChae et al.ll20l3 : 
ICappellari et ahl 20061). 


Emsellem et al. I 1 ^ 2011 ^ classify ATLAS®° ETGs by the 


degree of large-scale r otation using the criterion defined by 
lEmsellem et ahl (l2007l ). As noted in section 1, some previ¬ 
ous tests (based mostly on GC kinematics) of MOND in 
elliptical galaxies obtained contrasting results for fast and 
slow rotators [MOND was not successful for the outer part 
(> 2-3Re) kinem atics of slow rotators in several cases; see 
ISamuroW g (|2014| ) and references therein]. For SDSS ETGs 
this kinematic information is missing. The majority (86 ± 2 
per cent) of ATLAS®° ETGs are classified as fast rota¬ 
tors. However, for the above 11 ATLAS®^ ellipticals selected 
using essentially the same criteria as for the SDSS sam¬ 
ple, about one half (6 out of 11) are fast rotators. These 
two kinematic classes of galaxies have mean VDP slopes 
(rj) = —0.067 ± 0.014 (fast) versus —0.046 ± 0.017 (slow) 
where the discrepancy is not statistically significant based 
on the small samples. We caution, however, that the low 
precision of current samples does not allow a sensitive test. 

Fig. [2] shows the distributions of dark matter fraction 
within the effective radius (/dm) for the SDSS and the 
ATLAS^° nearly round galaxies as a function of stellar 
mass {Mi,), effective radius {Re) and Sersic i ndex (nser), 
where SDSS parameter values are taken from the lGhae et al.l 
||2014|) modelling result for their fiducial inputs. There are 
significant differences between SDSS and ATLAS^® samples. 
SDSS galaxies have larger /dm overall (mean (/dm) ~ 0.32 
versus 0.18) and are biased towards greater stellar masses 
(mean log]^g(M*/M 0 ) « 11.3 versus 10.8) and (more notice¬ 
ably) larger sizes (mean logiQ(i?e/kpc) « 0.82 v ersus 0.48). 
The r ange of modelling results considered by iGhae et al.l 
ll2014h taking into account various systematic effects give 
0-27 < /dm < 0.38. The second panel of Fig.[2]indicates that 
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Figure 3. Predicted distributions of VDP slope r] (equation 13) 
for ~ 2000 nearly round SDSS elliptical galaxies. Three cases are 
compared: Newtonian cases without and with DM presented in 
Chae et al. (2014) and a MOND case calculated here for (n, (/3)) = 
(1.5, 0) and shown in Fig. |4] Vertical red solid and dashed lines 
indicate the measured mean and its error for 11 ATLAS®® nearly 
round galaxies shown in Fig. [T] The arrows indicate the mean 
values of the distributions. 


the difference in (/dm) can be in large part attributed to the 
difference in galaxy sizes. Indeed, at log]^g(i7e/kpc) = 0.5 
SDSS galaxies have /dm ~ 0.25 (with a systematic error 
of ~ 0.05 due to input variations) getting closer to 0.18 of 
ATLAS®° galaxies. 

The difference in Re between the SDSS and the 
ATLAS®° samples may be due to various factors includ¬ 
ing different parent populations, measurement errors and 
different methods of modelling light profiles. The ATLAS®^ 
parent sample is a volume limited sample in the local Uni¬ 
verse wi thin D < 42 Mpc (alm ost) complete for M* > 6 x 
10 ®Mq (|GaDDeIlari et al1l201ll ). The SDSS parent sample is 
magnitude-limited to cosmological distances (up to redshift 
~ 0.2 with a median redshift of ~ 0 . 1 ) and thus completeness 
becomes less secure for less luminous galaxies and at larger 
distances. It is also likely that measurement/modelling er¬ 
rors are at work as can be seen, e.g., from the systematic dif¬ 
ference between RC3 an d 2MASS photometric data for the 
same ATLAS®^ galaxies l|GaDDelIari et al.ll201 j) and depen¬ 
dence on modelling details for SDSS galaxies i Meert et al.l 
l2013l) . 


We assume that the statistical distribution of rj is not 
significantly affected by such sample difference as that be¬ 
tween SDSS and ATLAS®° samples and use the measured 
distribution of rj for the ATLAS®® galaxies to interpret mod¬ 
elling results for the SDSS galaxies. We also make no dis¬ 
tinction between fast and slow rotators and note that our se¬ 
lected ATLAS®® galaxies are nearly evenly divided (Fig.[T]). 


3 RESULTS 

The output distribution of rj for the SDSS galaxy sam¬ 
ple can depend on the input stellar IMF (distribution), 
parameterised by 5m (equation 1121) . and VD anisotropy 
P{r) for a given MOND model (interpolating function). 
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V 

Figure 4. Predicted distribution of VDP slope r) (equation 13) for the SDSS galaxy sample depending on MOND interpolation function 
and VD anisotropy mean value (/3) assuming anisotropy is constant and has a galaxy-to-galaxy scatter of 0.2. We consider a MOND 
interpolation function of the form given by equation (5) with index n = 1, 1.5 and 2 as well as that given by equation (7) taken from 
Bekenstein’s (2004) relativistic theory of modified gravity. For each MOND model we consider three cases of (/3) = +0.2 (light gray), 
{/?) = 0 (blue) and (/3) = —0.2 (dark gray). Red lines are the same as in Fig. [3] The arrows indicate the mean values of the distributions. 


Our standard choice for th e VD-dependent stellar IMF is 
that preferred by ATLAS®° dCappellari et al1l2013lj) . which 
can be re-expressed in our adaptation as Sm = 0.11 + 
0.35logjQ(2“°'^®(Te2/130 km s“^). Unless specified other¬ 
wise, all presented results are based on the input of this 
ATLAS®° IMF distribution. However, as demonstrated be¬ 
low, it turns out that the output rj distribution is only weakly 
dependent on the input Sm distribution. 


We first co mpare the Newtonia n cases without and with 
DM presented in IChae et al.l (l2014l ) with a MOND case. The 
results for the three cases are shown in Fig. [3] The Newto¬ 
nian case with DM reproduced here is the result based on 
the fiducial inputs bv IChae et al.l (l2014l ). Newtonian dynam¬ 
ics clearly requires DM in a statistical sense (see also Fig. [3) 
and can match the obs erved VDP distri bution with adjust¬ 
ment of various inputs dChae et al.ll20lj ). However, it is also 
evident that MOND can reproduce well the VDP distribu¬ 
tion. 


We present in turn MOND results for constant 
anisotropy (section 3.1) and for varying anisotropy of the 
form /3dOM(r) (section 3.2). 


3.1 Constant anisotropy 


Anisotropy values are drawn from a Gaussian distribution 
with a mean of (/3) and a standard deviation of 0.2, which is 
motivated from studies of nearby early-typ e galaxies (e.g. 
Gnriiard et al.l I 2001 I : ICanpellari et al.] 20071 : iThomas et al.l 


200711 . The predicted distribution of r) for the SDSS galaxy 


sample then depends on the input mean anisotropy {j3) 
and MOND interpolating function. We consider (/3) = +0.2 
(light gray), 0 (blue), and —0.2 (dark gray) for each of n = 1, 
1.5, and 2 in the function given by equation as well as 
for the function given by equation 0. The results are dis¬ 
played in Fig. |4l For n = 1.5 the predicted mean of rj can 
match easily the measured value. The cases n = 1 and n = 2 
can also match the measured distribution but not for all the 
cases of (P). The case n = 1 prefers (/3) = +0.2 while n = 2 
prefers (P) = —0.2. These results imply that a number of 
degenerate sets of (n, (/?)) can be consistent with the mea¬ 
sured rj value. However, Bekenstein’s model is disfavoured 
for the assumed cases of constant anisotropies. 


In our approach it turns out that the output dis¬ 
tribution of rj is only weakly dependent on the input 
IMF. In the random process of finding {5m,P) where Sm 
(equation 1121) is essentially a representation of the stel¬ 
lar mass-to-light ratio depending on the IMF, the poste¬ 
rior distribution of Sm generally deviates from the prior 
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Figure 5. Input and output distributions of stellar IMF mis¬ 
match parameter 5m (equation 12). Three input IMFs are rep¬ 
resented by gray lines with an intrinsic scatter of 0.1 dex (rep¬ 
resented by dashed lines): Chabrier (2003), Salpeter (1955) and 
ATLAS^^ (Cappellari et ah 2013b). Data points are the output 
models and the black lines indicate the median relations (with 
dashed lines indicating the 68 percent limits). Thin lines rep¬ 
resent recently measured IMF behaviours: red - Conroy Sc van 
Dokkum (2012); green - Cappellari et ah (2013b); blue - Tortora 
et ah (2013). 


input distribution. Fig. [5] shows the distributions of 5m 
for three cases of in put IMF: two V D-independent cases of 
IChabrie^ (|2003l ) and Sal^ete d 19551) a nd one VD-dependent 
case of ATLAS^^ ( Cappellari et all 2013^), with an in- 
trinsic scatter of 0.1 dex ( Cappellari et all l2013bll im¬ 
posed in all cases. Irrespective of the input IMF the pos¬ 
terior mean of Sm increases with VD and its functional 
behaviour is well within the behaviours of recently in¬ 
ferred IMFs IConrov fc van Dokkumll2012l: Caja pellari et al.l 


l2013bl : iTortora. Romanowskv fc Napolitandl2013l l. 

For the cases of constant anisotropy the predicted width 
(containing 68 percent of galaxies) and standard deviation of 
rj (~ 0.23) are smaller than the measured standard deviation 
of arj « 0.037 ± 0.007. This discrepancy is largely removed 
if varying anisotropies are used as described in section 3.2. 



r/R, 

Figure 6. Examples of varying anisotropy y5dOM('^) = /^i/(l + 
-h /32/(1 + for mean (^e) = (/loo) = +0.2 (red), 0 

(black) and —0.2 (blue) and 0 < ri//2e < < 1- 


3.2 Varying anisotropy 

For ^dOM(r) = /3i/(l -I- rl/r"^) + /32/(l -f rijr^) the radial 
mean for 0 ^ r ^ he and the anisotropy at infinity 
Poo = Pi + P2 are assumed to be drawn from a Gaussian 
distribution with a standard deviation of 0.2. Because Pe 
can differ from Poo for a galaxy (although the means for the 
entire sample are assumed to be the same, i.e. (Pe) = (Poo)) 
and 0 < ri/Re < r^IRe < 1 the adopted function al¬ 
lows wild variation within Re. This does not necessarily 
mean that /3dOM(r) represents or fully encompasses real be¬ 
haviours of anisotropy. We just use /?dOM(r) to mimic some 
uncertainties in varying anisotropies. Examples (taken from 
modelling results) of /3dOM(r) for mean (/3e) = IPoo) ~ - 1 - 0 . 2 , 
0 and —0.2 can be found in Fig.|^ Similar examp les but with 
P(r = 0 ) 7 ^ 0 can be found in the appendix C of IChae et al.l 
( 120121 ) . 

Fig. [ 7 ] shows the predicted distributions of 77 for the 
same combinations of (n, (/3e)) as in Fig. [T] of the constant 
anisotropy case. Compared with the constant anisotropy 
case the predicted mean (77) shifts by ~ —0.01 for the same 
input of (n, (P)). Consequently, the empirical value of 77 is 
best matched by the case n = 1. The case tt = 2 is now 
disfavoured. Interestingly, Bekenstein’s model can now be 
consistent with the empirical value with a suitable choice of 
(Pe) although it may not be a preferred model. These results 
confirm that anisotropy shapes can matter in the study of 
MOND (as well as DM) in spheroidal galaxies. 

The width (containing 68 percent) of the distribution 
(~ 0.27) is about 20 percent larger compared with the con¬ 
stant anisotropy case, but still somewhat lower than the 
measured standard deviation. However, the marginal dis¬ 
crepancy may be a result of non-Gaussian nature of the dis¬ 
tribution. Indeed, the standard deviation of the distribution 
(~ 0.4) is larger than the 68 percent width and matches well 
the measured standard deviation. 
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Figure 7. Same as Fig. |4] except that radially varying anisotropies are used and /3e refers to the radial average of an anisotropy within 
Re of a galaxy. 


4 DISCUSSION AND CONCLUSIONS 


We have investigated VDPs for r < 7?e of ~ 2000 nearly 
round SDSS galax;ies with 1 0 < login(fU * /M(7)) < 12 (mean 
~ 11.3). Previous study by IChae et al.l ll2014h shows that 
the predicted distribution of VDP slope r) (equation 1131) 
cannot match an observed distribution without DM un¬ 
der Newt onian gravi t y (see also Fig. [3|. Furthermore, as 
shown in IChae et al.l (l2014l l each galaxy requires different 
amount and profile of DM. DM fraction within Re ranges 
from /dm ~ 0.1 - 0.5 (68% range: see Fig. ^ and the in¬ 
ner DM densi ty profile slope a, in the gener alised Navarro- 
Frenk-White jNavarro. Frenk fc Whitdl 19971 ') parameterisa- 
tion pDM(r) oc (r/rs)““(l-|-r/rs)“®+“, ranges from a « 0.7 
- 1.5 (68%). These results mean that elliptical galaxies do 
require DM under Newtonian gravity in a statistical sense, 
although some fraction of galaxies may not contain signifi¬ 
cant amount in the optical region. 

Our analysis of the same galaxies under MOND here 
indicates that the observationally-derived value of rj (in 
11 nearly round ATLAS®° galaxies) can be reproduced 
based on a single MOND interpolating function without 
DM applied to all the galaxies. It is particularly strik¬ 
ing that with an interpolation function similar to that of 
n = 1 or 2 in equation ®, used i n various MOND studies 
particularly for disc galaxies ( e.g. Fhmne^&_Binri^200^ 
Sanders fc Noordermeerl l2007l : lAngus et al.l l2008l : iMilgrori] 
2012 l ~). the empirical v value can be easily matched with a 
relaxed assumption of VD anisotropy (Figs. |4] and [7|). 


The specific functional form of MOND interpolat¬ 
ing function depends, however, on the assumption on VD 
anisotropy (Figs. [4] and [7|) and cannot be uniquely deter¬ 
mined from the current empirical rj value. The implied inter¬ 
polating function can be varied by varying mean anisotropy 
or radial behaviours. For constant anisotropy and the func¬ 
tional form given by equation mean anisotropy (/3) = 
+0.2, 0 and —0.2 give, respectively, 0.7 < n < 1.4, 1.1 < 
n < 1.9 and 1.3 < n < 2.5. For the case of radially varying 
anisotropy there is an overall shift of An ~ —0.5 for the same 
value of mean anisotropy so that the traditionally standard 
model n = 2 is disfavoured. These manifest a degeneracy 
between MOND interpolating function and VD anisotropy, 
which is reminiscent of the well-known mass-anisotropy de¬ 
gener acy under standard dynamics (e.g. ISinnev fc MamonI 
Il982ll . This interpolating function-anisotropy degeneracy 
can be alleviated by a more precise empirical value of {rf) 
or prior constraints on anisotropy. 

Our analysis has been limited to nearly round galaxies 
(under the spherical symmetry assumption) in the interme¬ 
diate acceleration regime although they cover a broad range 
of mass and size. A MOND analysis of ETGs of any elliptic- 
ity is considerably more challenging but would be necessary 
for a more rigorous test of MOND in spheroidal systems. 
If a random sample of elliptical galaxies (drawn from the 
same SDSS parent sample of ETGs using the same crite¬ 
ria except for the ellipticity limit) is analysed based on the 
MOND spherical Jeans equation (equation [^, the result- 
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ing » 7 -distribution has a similar mean but a somewhat larger 
scatter, probably hinting that MOND works also for ETGs 
of any ellipticity. 

Our analysis provides a novel and rigorous test of 
MOND in spheroidal galaxies in the sense that it is based 
on a statistically representative sample of spheroidal galax¬ 
ies and a range of possibilities of VD anisotropy. Our test is 
most sensitive at r < Re/2 where gravitational acceleration 
g due to the baryonic mass distribution is > uq (Fig. [8]). 
However, our galaxies cove r the range 0.1i?e ^ r < Re (see 
Fig. 1 of IChae et al.l 120141 ') for which O.Sao ^ 9 // lOOao 
(Fig. [8|. Our results support the view that MOND sug¬ 
gested by dynamics of rotating galaxies is also likely to be 
valid for dispersion-supported galaxies in their inner regions 
(r < Re). Our results are in line with the test of MOND 
through hydrostatic equilibrium in two X-ray bri ght ellip¬ 
tical galaxies over a wide acceleration range by iMilgromI 

ll2012ll . 

In this work we could not separately consider slow and 
fast rotators because the necessary kinematic information is 
not available for the analysed SDSS galaxies. The fact that 
the similarly selected ATLAS®® elliptical galaxies are nearly 
evenly divided (see section 2) hints that the SDSS sample 
may also contain both kinematic classes, although it is likely 
to be somewhat biased towards slow rotators as the galaxies 
are on average more massive and lar ger than th e ATL AS®® 
counterparts (see Fig. [2)). Recently. ISamurovid (|2014[ ) car¬ 
ried out individual Jeans modelling of ten elliptical galax¬ 
ies, four of whom are slow rotators, based on observed GC 
kinematic data up to several effective radii. Considering two 
cases of constant anisotropy /3 = 0, —0.2 and one case of 
radially varying /3(r) « 0.5r/(r -|- 1.4Re) for MOND inter¬ 
polating fun ctions of n = 1, 2 in equation m and that of 
equation amurovi3 (|2014| ~1 could not find successful fits 
of slow rotators without additional component of DM be¬ 
yond ~ 2-3Re. For one galaxy NGC 4486 (M87), which is 
in the central region of a cluster, the considered MOND 
models and anisotropies had difficulty of fitting GC kine- 
m atic data even in the inner region (r < 0.5i?e). However, 
as ISamurovi lim) notes, it remains unclear whether these 
problems based on limited cases of anisotropies imply break¬ 
down of MOND for certain objects. In this work we have 
considered a broader range of possibilities of anisotropy and 
found that VDPs of a mixed population of fast and slow ro¬ 
tators for r < Re could be statistically explained by MOND 
without any DM. 

If the empirical MOND law is to be truly meaningful, 
there must exist a universal interpolating function. Studies 
up to the present indicate that this is a reasonable possibil¬ 
ity. However, a firm conclusion can only be reached through 
precise determination of interpolating functions of galaxies 
of various types. Our analysis shows that the statistics of 
VDPs of elliptical galaxies provides a useful tool to test 
interpolating functions. This statistical analysis can com¬ 
plement individual analys es that can be pe rformed with de¬ 
tailed kinematic data (e.g. ISamuro^ lEHl). Despite its cur¬ 
rent limits, our analysis appears to support the simple in¬ 
terpolating function [equation ([5]) with n = 1] (see Figs. [4] 
and 0, which has also been preferred by galaxy rotation 
data (e.g. iFamaev fc Binn^l2005l : ISanders &: Noo rdermeerl 
l200if ) and studies of elliptical galaxies [see lMilgrom ( 2012l l 
and references therein], although we cannot yet rule out the 



logio(^/^e) 


Figure 8. Distribution of gravitational acceleration g{r) nor¬ 
malised by ao = 1.2 X s“^ due to stellar mass within 

radius r in ~ 2000 nearly round SDSS galaxies based on mod¬ 
elling results with n = 1, 1.5 and 2 for (/3) = 0 shown in Fig. |4] 
The solid curves represent the medians while the dotted curves 
represent the 95% limits. The observed VDs of SDSS galaxies 
analysed in this work are luminosity-weighted LOSVDs within 
projected radius Hapi which is peaked at Re/2 (represented by 
the vertical dashed line) but ranges from ~ 0.1i?e up to ~ Re- 


traditionally standard model with n = 2 or even the in¬ 
terpolating funct ion (eouation Pzll implied by Bekenstein ’s 
modified gravity llBekensteinll20M : IZhao fc Famaevll2006l ). 

If a single interpolating function turns out to explain dy¬ 
namics of all galaxies, that would have far-reaching implica¬ 
tions. In this respect it would be quite interesting to empiri¬ 
cally verify and determine the supposed universal MOND in¬ 
terpolating function with future data. For spheroidal galax¬ 
ies both statistical (for large samples, as was demonstrated 
in this study) and individual (for large radial extent, per¬ 
formed earlier and ongoing) analyses of VDPs will be useful. 

We are grateful to Srdjan Samurovic for providing a 
thorough and helpful report of the submitted manuscript. 
We also thank Maurice van Putten for useful discus¬ 
sions/comments and Mariangela Bernard! for helpful com¬ 
munication regarding SDSS galaxies. 
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